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Abstract 



Semiclassical transition probabilities characterize transfer of energy between "hard" 
and "soft" modes in various physical systems. We establish the boundary problem 



Ph: 

^ ' for singular euclidean solutions used to calculate such probabilities. Solutions are 



^ found numerically for a system of two interacting quartic oscillators. In the double- 

"X' 

;h ' well case, we find numerical evidence that certain regular minkowskian trajectories 



have approximate stopping points or, equivalently, are approximately periodic. This 
property leads to estimates of tunneling excitation probabilities in that system and 
suggests that similar estimates may be possible in other systems with tunneling. 



1 Introduction 

Very little is known at present about semiclassical transition probabilities in systems 
with more than one degree of freedom. Let x denote one of (possibly many) coordinate 
variables and \n) be the n-th energy eigenstate. Then, a typical semiclassical matrix 
element is {n2\x\ni) for \n2 — ni\ ^ 1. In what follows we always assume n2 > 
rii. This matrix element appears when one calculates the probability of exciting a 
system originally in state \ni) to state \n2) by applying a high-frequency external force 
proportional to x. Similar matrix elements appear in other problems where there is 
a transfer of energy between "hard" and "soft" modes. The states \n) then refer to 
the "soft" subsystem. In elementary particle physics, a typical problem of this sort 
is calculation of cross sections for production of many low-momentum particles in 
collisions of two high-momentum ones. Electroweak baryon number non-conservation 
at high energies is an example that was extensively discussed recently (for most recent 
reviews see ref.|]I|) and there are other interesting cases as well 0. 

In a system with one degree of freedom, semiclassical transition probabilities can 
be calculated by the method described in the Landau-Lifshitz textbook P|. A natural 
generalization of this method to many degrees of freedom is achieved by reinterpreting 
it in terms of singular classical solutions [Q . In ref . we have found some approximate 
singular classical solutions for self-interacting field theories in four dimensions. We 
have shown how they can be used to reproduce the factorial asymptotics of cross 
sections for multiparticle production at relatively low energies and how the crossover 
to a high-energy regime may occur. This approach was pursued further in ref. 0. 

There are several problems preventing a full-scale application of the method of 
singular solutions to field theory, (i) There is no formal derivation of the expression 
for transition probability used in @, Q; consequently, some interpretation problems 



arise as to exactly what quantity is being calculated, (ii) No exact singular solutions of 
the required form are known in systems with more than one degree of freedom (except 
when a reduction to one degree of freedom occurs) and it is not clear what are the 
existence and uniqueness conditions for such solutions, (iii) In general, one expects 
that exact singular solutions are accessible only numerically, and the appropriate 
numerical procedure needs to be formulated. 

The present work is an attempt to get at least partial resolution of these problems 
using the simplest non-trivial case of two interacting non-linear oscillators as an exam- 
ple. In Sect. 2 we describe a derivation of the double functional integral representation 
for the following transition probability 

W{E^,E^) = Y\{h\VE,xVE,\a)^ , (1) 

a,b 

where the sum is over a complete system of states and Ve is the projector onto the 
subspace of a given energy E. We then show that in the semiclassical approximation, 
the transition probability (^ is saturated by a singular solution on a certain contour 
in the complex time plane. We will comment specifically on the appearance of time 
dependence in the originally time-independent problem (|l]). In Sect. 3 we present 
the results of a numerical search for singular classical solutions in a system of two 
interacting quartic oscillators. We do not attempt a rigorous analytic resolution of 
the problem (ii) above but rather present numerical evidence for the existence and 
uniqueness of singular solutions under specified boundary conditions. The results 
of Sect. 3 refer to the case of a single- well potential. The double- well case, which is 
closer to the problem of electroweak baryon number violation, is considered in Sect. 4. 
We present there numerical evidence that the minkowskian solution that starts with 
zero velocity at a turning point of a periodic instanton (which is a periodic euclidean 
solution 1]^) has an approximate stopping point or, in other words, is approximately 



periodic. As we illustrate in Sect. 4, this property leads to straightforward numerical 
estimates of tunneling excitation probabilities in that double- well system. Similar es- 
timates may be possible in other systems with tunneling, maybe even in field theories. 
Concluding Sect. 5 contains a discussion of the results. 

2 Semiclassical formula for transition probability 

The projector Ve appearing in (^ has a simple representation, given in [^, only 
when the system reduces to a collection of non-interacting harmonic oscillators. In 
field theory, this usually happens at large positive and negative times when an initial 
configuration dissociates into free particles. A system of a few degrees of freedom, 
however, typically never gets out of the non-linear regime. To be able to use the 
representation of ref. for Ve in the case of two interacting oscillators, we need 
to approach the original problem in the following way. Two interacting oscillators 
represent a particle confined by a potential to some region on a two-dimensional plane. 
Let us couple this particle to electromagnetic field by giving it some small charge. 
Then, an excited state of the particle will eventually decay into photons, for which 
our formula for the projector is applicable. For a highly excited state, this decay can 
be viewed as an almost classical electromagnetic radiation by the oscillating particle. 
So, the coupling to electromagnetic field does not spoil the semiclassical nature of our 
problem. Moreover, because the radiation takes place in real time, the action on the 
corresponding part of the solution does not contribute to the exponential factor in 
the probability (|I]). Thus, making the electromagnetic coupling arbitrarily small, we 
make the probability arbitrarily close to that in the system without electromagnetic 
coupling, even though at large times the two systems look entirely different. 

With this modification in mind, we can use the formula for Ve at sufficiently large 



times. On the other hand, the states in eq.(|I]) are taken at some fixed moment of 
time t = to- To proceed, we should define the state VEila) at t = — oo and the state 
'PE2W) at t = 00, where we know how to do that, then evolve these states to t = to 
and use the result in eq.(|lD. As usual, it is convenient to define initial and final states 
in the interaction representation. Then, eq.(|l]) takes the form 



W^(Ei,E2) =E (Koo)|Pi^,^(oo,to)x^(to,-oo)PsJa(-oo)) , (2) 



b 
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where 5'(cx2,to) is the "half of the S'-matrix" operating from t = to to t = 00. Thus, 
in order to use the known expression for the projector we introduce time dependence 
in the originally time-independent problem. This may seem to be an additional 
complication but in fact it is just the opposite. Instead of calculating difficult overlap 
integrals of wave functions, we find the exponential suppression that is present in 
W{Ei, E2) simply from the imaginary part of the action of a classical solution. 

In the leading semiclassical approximation, we are interested only in the exponen- 
tial factor in the probability, so the operator x in eq.(0) is inessential. (Of course, 
without this operator the pre- exponent will vanish.) The exponents of the two halves 
of the S'-matrix in @) add up to that of the whole S-matrix. We can now proceed in 
a rather close analogy with the calculation of ref. 0]. The main difference with that 
calculation is in the contour in the complex time plane, on which a classical solution is 
to be found. When we consider a transition probability rather than a transition am- 
plitude, we actually need a contour that consists of two parts symmetric with respect 
to the real axis. In the present case, we use the contour shown in Fig.l together with 
its reflection into the lower half-plane. Unlike the contour used in 0, the contour of 
Fig.l has both euclidean and anti-euclidean parts, in addition to minkowskian ones. 
A classical solution along the contour shown in Fig.l describes a transition from a 
state of energy Ei ai t = — cxd + iTi/2 to a state of energy E2 a.t t = +cxd + iT2/2. 
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Figure 1: Contour in the complex time plane, on which we look for a real classical 
solution. 

This solution is necessarily singular because the initial and final states have different 

energies. By convention, the singularity is positioned at t = 0. 

It is important to realize how crucial is the fact that we consider a total probability 

and therefore sum over \a) and \b) in eq.(^. If we were to consider probability of 

transition between certain \a) and \b), we would have to impose the corresponding 

boundary conditions on the real axis, at t = — cxd and t = +00. Shifting the points 

at infinity into the complex plane would then require an adjustment in the boundary 

values of the fields, leading to complex solutions. A complete system of states, on the 

other hand, can be inserted anywhere in the complex time plane, so by considering 

the total probabihty (g) we allow ourselves to restrict attention to real solutions]] 

* The question remains, which we did not attempt to resolve in this paper, whether additional 
singularities that a solution can have in the complex time plane are crossed when we deform the 
contour into the shape shown in Fig.l. One can verify, however, that no such singularities have to 



The derivation of the functional-integral representation for the probability (0) 
using the representations for the S-matrix and the projector Ve given in is now 
completely straightforward and we omit the lengthy intermediate expressions. The 
leading semiclassical result for the transition probability has the form 

W{Eu E2) ~ exp (EiTi - E2T2 - 2S) = exp(-Q) , (3) 

where S is the net euclidean action acquired on the euclidean and anti-euclidean parts 
of the contour of Fig.l. The action S corresponds to a solution that at t = iTi/2 
starts from a turning point with energy Ei, follows euclidean evolution to a singularity 
at t = 0, then follows anti-euclidean evolution to a turning point with energy E2 at 
t = iT2/2. The actions corresponding to the euclidean and anti-euclidean evolutions 
separately are infinite but they come with opposite signs and when added together 
give a finite net amount S. 

The minkowskian parts of the contour do not contribute to eq.(|^). Therefore, in 
applications of eq.(^ we can forget we ever needed to couple our system to the photon 
"bath" because it is always possible to choose the coupling to photons so small that 
it does not influence the euclidean and anti-euclidean parts of the evolution. 

As we will see in the next section, there are infinitely many singular solutions of 

the form described above, corresponding to infinitely many possibilities to choose the 

initial turning point at t = iTi/2. However, there is a preferred one, for which Q in 

(^ attains the smallest possible value at given Ei and E2. This solution saturates 

the semiclassical transition probability and determines the states at t = —00 -|- iT\j2 

and t = -|-oo -|- iT2/2. It also determines Ti and T2 for given Ei and E2. 

An analogous calculation in field theory will have to impose a momentum cutoff 
be crossed in the case of a single symmetric quartic oscillator where singular solutions can be easily 
found analytically. 



that specifies what modes are considered "soft". This is not an ambiguity of the 
method because such a cutoff has to be buih in any procedure of identifying "soft" 
particle production in experiment. One might think then that the transition proba- 
bihty is saturated by a solution for which the final state \b) has particles only in a 
few "hardest" of the remaining "soft" modes, because this allows for smaller multi- 
plicity. However, it is not so clear what the actual answer is: the tendency towards 
smaller multiplicity has to compete with the tendency to have a spatially localized 
coherent field configuration, which enhances the probability by achieving stronger 
non-linearity. Our computations for two oscillators described in Sect. 3 show that 
the tendency towards stronger non-linearity is indeed rather effective in distributing 
energy between "harder" and "softer" modes in the final state. 

3 Numerical results 

We now turn to numerical results for the case of two interacting quartic oscillators. 
The hamiltonian is 

H = ^x' + ^f + Vix,y), (4) 

where the potential V{x,y) is 

V{x, y) = ]^x^ + hoY + \{^' + 6x2y2 + y^^^ . (5) 

The particular choice of the ratios of couplings in (^ is motivated by the property 
that the system (^-(§) is equivalent to the A0^ theory on a space "lattice" consisting 
of two sites. The euclidean and anti-euclidean parts of classical evolution can be 
viewed as mechanical motions in the potential equal to —V{x^y). 

Instead of the coordinates x and y it is convenient to use two other variables - 
the polar angle (f) in the (x, y) plane and the value of the potential V . The euclidean 



part of a singular solution starts with zero velocity at some angle (pi and V = Ei ^ 0, 
and the anti-euclidean part stops at angle 0/ and V = E2 > Ei. The main question 
is how we match these two parts at the singularity. 

To formulate the matching condition, consider the following limiting procedure. 
Let us choose some large value of the potential A ^ E'l, £'2- The euclidean trajectory 
starting with zero velocity a.t(j) = (f)i,V = Ei will reach, in some time T'i(A), the value 
y = A at some angle = 0a- The euclidean action accumulated on this trajectory 
is some 5*1 (A). We now look for the trajectory that starts with zero velocity at some 
angle 0/(A) and V = E2 and reaches, in some time T2(A), exactly the same point 
V^ = A, = 0A as the first one. Mechanical intuition suggests that such trajectory 
always exists and our numerical results also show so. The corresponding euclidean 
action is some 5*2 (A). Due to anharmonicity, both trajectories reach infinite values of 
V in finite times, in other words, limits 

Ti,2 = lim Ti.2(A) (6) 

always exist. A far more subtle property is the existence of limits 

S = hm (5i(A)- 52(A)) , (7) 

A— >oo 

0; = l^M^)- (8) 

If these limits exist, they define the singular solution corresponding to given Ei, E2 
and a given initial configuration 0j. We then scan among all initial configurations of 
energy Ei and find the one that leads to the smallest suppression factor in (^. To 
summarize, the euclidean and anti-euclidean parts of a solution are matched at some 
large value of the potential energy A and then limit A — > 00 is taken in the hope that 
the quantities d^)-® remain stable. In what follows, we present numerical evidence 
in favor of the existence of the limits (0)-(§). 

9 



Numerically, matching at finite A is simplified by the following property of eu- 
clidean trajectories, which we found empirically. Suppose there are three trajectories 
starting with V = E and zero velocity at angles 0i, 02 and 03 and passing through 
y = A at angles $1, $2 and $3, respectively. The property is that if 0i < 02 < 03, 
then $1 < $2 < *^3- This "monotony" allows us to use the familiar root finding 
methods of bisections and secants, to find the trajectory that starts at V = E2 and 
passes V^ = A at a specified value of the angle. 

For our computations we chose the following values of the parameters: u = 1.5, 
A = 0.1, El = 0.0001. Such a small value of Ei means that the system is excited 
out of its ground state, which is the case of main interest in many applications. We 
consider a set of different values of the final state energy E2 in order to obtain the 
dependence of the suppression factor @ on the excitation energy. 

In a numerical computation, it is of course impossible to make A arbitrarily large. 
Let us estimate the error introduced by the finiteness of A. In the Landau-Lifshitz 
formula for a quartic oscillator, with V{x) = x^/2 + gx^/A, the absolute value of the 
error in Q/2 introduced by the finiteness of A is 

v^£ {(Vix) - Elf- - (Vix) - E^fl^) dx = ^^ + 0(A-3/4) (9) 

Guided by eq.(^, we estimate the same error for our system (assuming the limits 

(0)-(D exist) as 

E'} — £/i Eo — E^ 

<-|7TT7^, (10) 



(A,ffA)V4 - (^A)i/4 

where Agg is an effective coupling constant for a given region of angles, and the 
inequality follows from the fact that the smallest anharmonicity is along the x = 
and y = axes, and this corresponds to Agg = A. 



If we want the error ( |10D to be smaller than some e, it is sufficient to take A equal 
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to or larger than 



A.4(^y. 



Large values of A dictate the need for high accuracy in solving the euclidean evolution 
equations and matching euclidean and anti-euclidean parts of a solution, in order to 
insure that numerical errors coming from those sources do not add up to quantities 
larger than e. Suppose the numerical error we make in solutions for x{t) and y{t) is 
some 6. The abbreviated actions 5*1,2 (A) — -£'i,27'i,2(A)/2 of the euclidean and anti- 
euclidean parts of a solution are found according to 

/■Ti,2{A)/2 

5i,2(A) = 2 / {V[x,,2{t),yi,2{t)] - Ei,2}dt . (12) 

Jo 

Apart from solutions staying on the x = and y = axes, for which our system 

reduces to that of one degree of freedom, every euclidean solution eventually reaches 

a region where the absolute values of both x and y are large. This is the region where 

the main error in Si, 2 (A) comes from. In such a region, the system (^ with our 

particular choice of the potential, eq.(^, is effectively separable in terms of variables 

^ = (x + y)/^/2 and r] = {x — y)/\^, so the main error in S'i,2(A) is the sum of errors 

in corresponding abbreviated actions. The errors in ^ and r] are not larger than \/2S. 

Then, for example, the error in the potential for ^, V^, is \/25dV^/d$, and the error in 

the abbreviated action for ^ is estimated as 



2V25J ^dt-2Sj jy^^^^ 



= 46^V^{U) + 0{6/^V^{U)). (13) 

(Strictly speaking, some numerical error comes also from determination of Ti 2 but 
that can be easily made much smaller than (0).) Even though both V^{C,a) and 
Vri{ri\) are large, only one of them is of order (and hence essentially equal to) A. So, 
the numerical error in each of the abbreviated actions ([T^), produced by the error 6 
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A 


0.5 X A, 


A, 


2.5 X Ae 


5 X A, 


10 X Ae 


g/2 


23.05 


23.06 


23.08 


23.09 


23.10 


H 


1.40979841 


1.40979916 


1.40979996 


1.40980045 


1.40980085 



Table 1: Values of half of the exponent Q of eq.(3) and the final state angle 0/ for 
E2 = 15, initial state angle (pi = 7r/6 and different values of the matching energy A. 
A, = 5.06 X 10^ 



in the solution, is estimated as A6yA. If we want the corresponding error in Q/2 to 
be much smaller than e, we have to make S much smaller than 

5, (A) = -^ . (14) 

For the smallest admissible A, eq.(0), this becomes 

We see that even moderate improvements in the desired accuracy e of Q/2 require 
considerable improvements in the accuracy of the solution. The required accuracy of 
the solution is also higher for transitions between states that are more separated in 
energy. 

Our computations were done with e = 0.1. The equations of motion were solved 
and euclidean and anti-euclidean parts were matched to accuracy more than an order 
of magnitude better than (|1^. Table 1 illustrates the existence of the limits (l^)-® 
for E2 = 15 and (pi = it/6. Note that for A > A^, the value of Q/2 does not change 
with A within the accuracy e, in correspondence with our estimates. We remind that 
each value of Q/2 in Table 1 is obtained as a difference between large (divergent at 
A -^ 00) abbreviated actions Si — EiTi/2 and 5*2 — E2T2/2 acquired on the euclidean 
and anti-euclidean parts of the solution. For example, for A = 5 x V^, the values of 
these abbreviated actions are of order 6 x 10^. 

12 



4>i 


~ 10-19 0.008 


0.017 0.028 


0.041 


0.058 


0.083 


0.129 


0.231 


0.537 


0/ 


7r/20 


vr/lO 37r/20 


7r/5 


7r/4 


37r/10 


77r/20 


27r/5 


97r/20 


Q/2 


25.2 23.0 


21.9 21.3 


21.1 


21.0 


21.2 


21.6 


22.2 


23.1 



Table 2: Several values of the initial angle and the corresponding values of the final 
angle and Q/2 for E2 = 15. Estimated error in Q/2 is less than 0.1. 



E2 


10 


15 


25 


50 


Q/2 


15.4 


21.0 


30.9 


51.8 



Table 3: Values of Q/2 determining suppression factors (3) for initial state energy 
El = 0.0001 and various final state energies E2. Estimated error in Q/2 is less than 
0.1. 

Next we consider values of the exponent Q for fixed E2 and different values of 
the initial angle 0j. Results for E2 = 15 are presented in Table 2.[| We see that Q 
has a minimum near (pf = tt/A. The corresponding value Q/2 = 21.0 determines the 
suppression factor (H) for E2 = 15. 

Note that the preferred final state is not purely the "harder" ?/-component but 
a state where x and y are comparable, which achieves stronger non-linearity. This 
tendency persists to larger values of 00. For example, aX. uj = 3, E2 = 10, very few 
?/-quanta need to be produced but the system prefers 0/ near 37r/10 rather than 
<Pf = 7r/2. 

The minimum in Q/2 becomes more pronounced as the final state energy E2 
increases. Table 3 gives such minimal values of Q/2 for various values of i?2- The 

suppression increases with the excitation energy, as may be intuitively expected. 

t These results were actually obtained by finding (f>i for a set of (fif equally spaced between and 
7r/2. Equal spacing of (j)f, rather than of </>,;, happens to give more points in the most interesting 
region of angles where Q is close to its minimum. 
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4 Tunneling probabilities 

In the previous sections, we considered transition probabilities, for which we summed 
over all initial and final states of given energies (or, more precisely, states with energies 
distributed in a small interval near given ones - microcanonical states). For potentials 
with two or more local minima, at energies less than the height of the barrier, we 
may be interested separately in the probability of a tunneling process, when a system 
originally distributed near energy Ei in one well gets excited to the state near energy 
E2 in another. If Ei = E2, tunneling probabilities for microcanonical initial states 
can be calculated using regular euclidean solutions - periodic instantons 0. When 
El ^ E2, however, we again have to use singular solutions. 

In the case of a single degree of freedom, the relevant singular solution is easily 
found (see for example ref.pj). Let us assume Ei < E2. The periodic instanton 
corresponding to energy E2 has two turning points X2 and X3, one of them, X2, in the 
well where the system is originally located. If we use X2, with zero velocity, as an initial 
condition for minkowskian evolution, the resulting (real) minkowskian trajectory will 
stop again at a different point xi in the same well. This new turning point xi, which 
still has energy E2, can now be used as the final point of the anti-euclidean segment of 
a singular euclidean solution whose euclidean segment starts at some xo with energy 
El. The whole trajectory from xq to 2:3 comprises the singular solution that describes 
the tunneling process. The corresponding suppression exponent is the sum of the 
exponent Q, eq.(||), acquired from xq to xi, and the suppression exponent associated 
with the periodic instanton. 

The stopping of the minkowskian trajectory at Xi is a special feature of systems 
with one degree of freedom. When there are several degrees of freedom, velocity has 
several components which in general will never vanish simultaneously. Further, there 
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seems to be no reason why they should vanish simultaneously (except at the initial 
point) for the minkowskian trajectory starting at the turning point of a periodic 
instanton. However, a numerical "experiment" shows that in a double- well system of 
two interacting quartic oscillators (see below) this in fact nearly happens. In other 
words, the minkowskian solution starting at the turning point of a periodic instanton 
is approximately periodic. 

Because the periodicity of the minkowskian solution is not exact, the method de- 
scribed above for the case of one degree of freedom cannot be literally applied to 
several degrees of freedom. Most likely, to generalize this method to several degrees 
of freedom, one has to consider complex singular solutions. The approximate pe- 
riodicity is still of value, however, because it allows us to construct a purely real 
approximate singular solution describing a tunneling transition and use it to estimate 
the probability of that transition. 

Specifically, consider the hamiltonian (|^) with the following double-well potential, 

V{x, y) = -^x^ + ^uV + ^{^' + Q^V + y') - oiy . (16) 

The term —ay was added to avoid the trivial situation when stopping occurs simply 
because y = on the whole trajectory. Below we present results for tu = 1.5, A = 0.01, 
a = 1 and E2 = —15 (minima of the potential (|1^) have negative energy). Similar 
results were obtained for other values of the parameters and for other quartic double- 
well potentials. The approximate stopping property is less pronounced for potentials 
of higher order. 

We found it convenient to characterize configurations in the left well by their polar 
angles (p in the polar system in which the origin is placed at a; = — 1/A^/^, y = 0. (This 
would be the location of the left minimum of the potential at a = 0.) The turning 
point of the periodic instanton for the quoted values of the parameters is = 0.0817. 

15 



Figure 2: (a) Absolute value of velocity on the minkowskian trajectory starting 
with zero velocity at a turning point of a periodic instanton. (b) The same for an 
arbitrarily chosen initial point. 

The absolute value of velocity v = (x^ + y^)^/^ on the minkowskian trajectory starting 
at this point with zero velocity is plotted as a function of time in Fig. 2a. We see that 
the absolute value of the velocity almost reaches zero. Numerically, the minimum 
absolute value of the velocity (not counting zero value at t = 0) in this case is about 
3 X 10^^, which is at least hundred times smaller than its "natural" value. For 
comparison, Fig. 2b shows the same quantity for the minkowskian trajectory starting 
with zero velocity at arbitrarily chosen point = 0.01. Here the absolute value of 
velocity does not approach zero (except at zero time) and is clearly aperiodic. 

We can use the approximate stopping point of a minkowskian solution as the 
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E2 


-24 


-15 


-10 -5 -1 


-0.25 


Qtun/2 


87.6 


68.5 


59.4 50.8 44.3 


43.0 



Table 4: Values of Qtun/^ determining estimated suppression factors (10) for the 
potential (9), initial state energy Ei = —24 and various final state energies E2. 
Estimated numerical error in <5tun/2 is l^ss than 0.1. 

final state 0/ for a real singular euclidean solution similar to those found in Sect. 3. 
Because the euclidean solution reaches exactly zero velocity at this point while the 
minkowskian solution does so only approximately, the trajectory comprised of the two 
will not be an exact solution. Still, it is useful for estimates of tunneling probabilities 
in our system. We thus describe the tunneling process as consisting of two stages: 
first, the system is excited from energy Ei to energy E2 while staying in the same 
well, then it tunnels at energy E2. The tunneling probability is then estimated as 

iytun(^i, E2) ~ exp(-Q + E2T' - S') ^ exp(-gtun) , (17) 

where Q is associated with the singular euclidean solution in the same way as in 
eq.(|^), S' is the action of the periodic instanton and T' is its period. Table 4 shows 
values of Qtun/2 for the potential (|T^) with the same values of u, A and a as above. 
El = —24 and various values of E2. For reference, the energy of the saddle point 
separating two minima is Eq = —0.222. We see that the suppression decreases as this 
energy is approached, though it does not go away completely. 

5 Conclusion 

Our aim in this work was to introduce the boundary problem for singular euclidean 
solutions and to gain evidence for its correctness by doing numerical calculations for 
a system of two interacting non-linear oscillators. For two oscillators with a quartic 
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double-well potential, we also presented numerical evidence that the minkowskian tra- 
jectories starting with zero velocity from turning points of periodic instantons (which 
are periodic euclidean solutions) have approximate stopping points or, equivalently, 
are approximately periodic. This property allowed us to construct real approximate 
singular solutions describing tunneling transitions and use them to estimate the prob- 
abilities of such transitions. 

One may try to extend the latter approach to other systems, including field the- 
ories. The accuracy to which the approximate stopping property holds will differ for 
different systems, and so will the accuracy of estimates based on it. However, within 
that accuracy, this approach gives an in principle straightforward method, essentially 
identical to the method used for a single degree of freedom, of numerical estimates of 
tunneling excitation rates in systems with many degrees of freedom. 

The author is grateful to V. Petrov, V. Rubakov and P. Tinyakov for useful dis- 
cussions, and to the Aspen Center for Physics, where some of these discussions took 
place, for hospitality. Part of this work was done while the author was at the Uni- 
versity of California, Los Angeles as a Julian Schwinger postdoctoral fellow. 
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